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ABSTRACT 

We argue that jittering jets, i.e., jets that have their launching direction 
rapidly change, launched by the newly formed neutron star in a core collapse 
supernova can explode the star. We show that under a wide range of parameters 
the fast narrow jets deposit their energy inside the star via shock waves, and 
form two hot bubbles, that eventually merge, accelerate the rest of the star 
and lead to the explosion. To prevent the jets from penetrating through the 
collapsing stellar core and escape with their energy, instead of forming the hot 
bubbles, the jets should be prevented from drilling a hole through the star. This 
condition can be met if the jets' axis rapidly changes its direction. This process 
of depositing jets' energy into the ambient medium is termed the penetrating jet 
feedback mechanism. The feedback exists in that the neutron star (or a black 
hole) at the center of the core collapse supernova shuts off its own growth by 
exploding the star. The jets deposit their energy at a distance of ~ 1000 km from 
the center and expel the mass above that radius. In our model, the material near 
the stalled shock at several hundreds kilometers from the center is not expelled, 
but it is rather accreted and feed the accretion disk that blows the jets. The 
neutrinos might influence the accretion flow, but in the proposed model their 
role in exploding the star is small. 



INTRODUCTION 



Jets can play a key role in exploding core collapse supernovae (CCSNe; e.g., Khokhlov 
et al. 1999; MacFadyen et al. 2001; Maeda & Nomoto 2003; Woosley & Janka 2005; Couch 
et al. 2009). Khokhlov et al. (1999) injected a jet with a radius of r^o = 1200 km at a 
distance of -Rin = 3820 km from the center. The mass in the two jets was ~ O.IM©, and 
their speed ~ 0.1c. Practically they injected slow massive wide (SMW) jets to explode the 
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star. As the jets launched by the newly formed neutron star (NS) are hkely to be narrow 
and fast (a velocity of > 0.1c), the formation of the jets launched by Khokhlov et al. (1999) 
should be explained. 

Couch et al. (2009) found that to match observations their jets must start with a large 
fraction of thermal energy (their model vlml2). The total mass in the jets in their model 
vlml2 was O.I2M0, and the maximum velocity 10^ km s~^. As evident from their fig. 7, 
such jets practically start as wide jets. Namely, their initial conditions for their model vlml2 
was that of a SMW outflow. Again, the formation of the SMW outflow requires explanation. 

MacFadyen et al. (2001) injected jets at a radius of Ri^ = 50 km, but their jets were 
injected at a much later time in the explosion, and are less relevant to the present paper. 
In any case, they also showed that SMW jets are efficient in removing the envelope further 
out. When narrow jets are simulated (e.g.. Alloy et al. 2000; Zhang et al. 2003, 2004) 
no envelope ejection occurs. Indeed, narrow jets that maintain a constant direction do not 
expect to expel a non-turbulent surrounding gas (Sternberg et al. 2007). In another set of 
simulations, Zhang et al. (2006) added a SMW outflow to their simulations of relativistic 
jets. The papers cited above further emphasize the need to answer the question as of how 
such SMW jets that can explode CCSNe are formed. 

Soker (2010; hereafter paper 1) argued that SMW jets in CCSNe, as simulated by, e.g., 
Khokhlov et al. (1999) and Couch et al. (2009), can be formed when fast {vj > 0.3c) and 
hght (total mass of ~ 0.01 — O.OSMq) jets interact with the infalling gas at distances 
of (~ 10^ < r < 10^ km) from the center. In paper 1 the same mechanism that was used to 
explain the formation of SMW jets in other astrophysical systems was proposed to operate 
in CCSNe. This mechanism is the penetrating jet feedback mechanism, that was applied to 
the growth of the super massive black hole during galaxy formation (Soker 2009; Soker & 
Meiron 2010), and to stellar environments (Soker 2008). When the jets penetrate through the 
ambient medium to large distances, they don't heat and don't expel the gas, and accretion 
onto the central region might continue. When the jets don't penetrate, but rather deposit 
their energy in the inner regions, they heat and expel the gas, assuming they are energetic 
enough and the hot bubbles don't suffer substantial non-adiabatic cooling. 

Maeda & Nomoto (2003) studied hypernovea explosions and injected bipolar jets at a 
radius of rjo — 1000 km for 1 sec, inside either a 25Mq or a AOMq star. The jets were 
injected at a constant direction. Their jets had a half opening angle of either 15° or 45° and 
a kinetic energy in the range of 1 — 30 x 10^^ erg. Their results seem to indicate that the 
stellar mass near the equatorial plane (perpendicular to the jets' direction) is not expelled 
efficiently, and the observed explosion is expected to be highly non-spherical. The jets we 
assume to be launched by the newly formed neutron star are similar in properties to the 
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jets simulated by Maeda & Nomoto (2003). However, there is a substantial difference. We 
assume the jets' axis to constantly vary, such that the envelope is expelled in all directions. 

For the efficient operation of the penetrating jet feedback mechanism in expelling the 
ambient gas, the relative direction between the narrow jets and the ambient gas must be 
changing, such that the jets do not manage to drill a hole and escape with their energy. 
The jets are rather shocked inside the ambient gas and heat it and/or expel it. The change 
in direction between the jets and the ambient gas is achieved either by transverse motion 
of the ambient gas, or changes in the jets' launching direction, e.g., an ordered precession 
(Sternberg & Soker 2008) or random direction changes termed jittering. Yet another way to 
prevent penetration is by wide jets (Sternberg et al. 2007). In that respect, the formation of 
wide jets by the newly formed NS or black hole, as found by Dessart et al. (2008), favors the 
operation of the model proposed here. In the magnetic field model of Dessart et al. (2008) 
the angular momentum is large, and the jets are expected to maintain a constant direction, 
and might have difficulties to expel the gas perpendicular to their direction. Adding a little 
jittering (see section |3]) will enable the jets to explode the entire envelope. The model studied 
by Dessart et al. (2008) requires high initial angular momentum. In the model discussed 
here there is no such requirement. 

Winds blown by the accretion disk around the newly formed NS have been also consid- 
ered in the literature. MacFadyen & Woosley (1999) found a disk-wind in their simulations, 
that has the same properties as SMW jets. However, the disk wind blown from an extended 
disk surface is less efficient (energetically speaking) than SMW jets formed from shocked 
narrow jets that are launched from the very inner region of the accretion disk. The reason 
is that a wind from an extended disk region is launched from a shallow potential well. In 
the penetrating jet feedback mechanism the explosion is powered by narrow jets that are 
launched from the inner disk, r ^ 30 km, where the potential well is much deeper than at 
the rest of the disk, and the accreted gas releases much more gravitational energy. Indeed, 
MacFadyen & Woosley (1999) simulated the formation of a black hole (and not a NS), and 
the strong disk wind is formed only when the collapsing core is rapidly rotating. 

Kohri et al. (2005) conducted a study of disk wind in CCSNe, where the central object 
is a NS. They propose that the wind energy is able to revive a stalled shock and help to 
produce a successful supernova explosion. Here again, the wind comes from an extended 
region in the disk, and it is less efficient than fast jets blown from the very inner region 
of the disk. Moreover, to form their proposed disk wind, Kohri et al. (2005) required the 
progenitor's core to rotate very rapidly, as they form the accretion disk earlier than in the 
present model. Another difference is that in our proposed model the jets deposit their energy 
further out, and we let the material near the stalled shock to be accreted. This is a major 
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difference between our model and models that are based on neutrinos. Namely, that our 
model does not revive the stalled shock. To the contrary, we need the material there to be 
accreted and form the accretion disk that launches the jets. 

In paper 1 a transverse motion of the infalling gas was considered, and a simple cooling 
by neutrino formula was applied. In the present paper we consider jittering jets, and show 
that for a wide range of parameters these jets can deposit their energy at r ~ 10^ km and 
form hot bubbles. These high-pressure bubbles accelerate the rest of the star and explode it. 
Before the bubbles merge they accelerate the material in their vicinity and form the SMW 
outflow. As the proposed mechanism is different in some key issues from previous models 
that are based on jets and winds from the NS vicinity, we present the basic ingredients of 
the model in sections El and El In section HI we use the penetrating condition to find the 
radius where the jets deposit most of their energy. Neutrino cooling is discussed in section 
ini and our summary is in section [6l 

2. THE PENETRATING JET FEEDBACK MECHANISM 

The basic ingredients and assumptions of the proposed explosion mechanism were de- 
scribed in paper 1. Here we concentrate on jittering jets, and use the following, somewhat 
different, assumptions and processes. 

1. We assume that about a Chandrasekhar mass of the iron core ~ 1.2Mq has already 
collapsed, and most of it already formed the almost final neutron star (NS). Namely, 
a mass of ~ I.IM© and a radius of ~ 20 — 30 km is accreting the rest of the inflowing 
mass. The corresponding escape velocity is fesc — 10^ km s~^. 

2. Based on recent numerical simulations (Ott et al. 2009; Blondin & Mezzacappa 2007; 
Nordhaus et al. 2010; Rantsiou et al. 2010) we assume that the angular momentum 
of the accreted gas onto the NS changes stochastically, and with a large enough ampli- 
tude to form an intermittent accretion disk with a rapidly varying angular momentum 
direction. Based on their simulations, we find that during a large fraction of the time 
the specific angular momentum allows for the formation of an accretion disk at a radius 
of ~ 30 km, and that the angular momentum direction is changing on a typical rate 
of ~ 10 rad s^^. In addition to the stochastic behavior behind the stalled shock at 
~ 300 km (Ott et al. 2009; Blondin & Mezzacappa 2007; Rantsiou et al. 2010), the 
magnetic axes of many NS are known to be inclined to the spin. Over all, the jets' 
direction is not expected to be constant, unless the core of the collapsing star has a 
large amount of angular momentum. 
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3. The properties of jets launched by NS (or BH) have some universal properties, such as 
the fast jets' speed is Vf ^ Vcsc — 10^ km s~^. 

4. There is a universal ratio between mass lose rate in the two jets to mass accretion 
rate. Using the same ratio of ejected to accreted energy as in Soker (2009), 0.05, 
the ejected (into both jets combined) to accreted mass average ratio is taken to be 
7] = Mj/Macc — 0.1. This average ratio includes the time periods when there is no 
accretion disk, or the disk does not blow jets. 

5. The ambient density profile in the specific calculations presented here are taken to be 
a power law, based on existing calculations (Wilson et al. 1986; Mikami et al. 2008) 
of a 15Mq core collapse star, rather than from the assumption of an inflow used by 
Soker (2010). The mass inside that radius of r < 300 — 1500 km that is accreted and 
form the accretion disk after bounce amounts to ~ 0.1 — O.5M0. 

From these assumptions we derive the typical kinetic energy carried by the jets 



With an efficiency of kb ~ 50% to channel the jets' energy to explosion, these jets can account 
for the kinetic energy of exploding CCSNe. 

If the jets penetrate through the surrounding gas they will be collimated by that gas, 
and two narrow collimated fast jets will be formed, similar to the flow structure in the 
simulations of Sutherland & Bicknell (2007) for AGN jets, and of MacFadyen et al. (2001). 
If, on the other hand, the jets cannot penetrate the surrounding gas they will deposit their 
energy in the inner region. A pair of two hot bubbles is formed. The bubbles will merge 
with each other and with earlier and later pairs of bubbles. Each bubble will accelerate the 
material around it over a very large solid angle. The wide outflow has more mass than the 
originally narrow jets, and by energy conservation the wide outflow has a lower velocity. 
A slow massive wide (SMW) outflow (jets) has been formed. It is this SMW outflow that 
expels the rest of the stellar gas. That a wide outflow can expel the stellar gas has been 
studied analytically by Kohri et al. (2005), and was demonstrated in numerical simulations 
by Maeda & Nomoto (2003) and Couch et al. (2009; where their high thermal energy jets 
are practically SMW jets). The results here offer an explanation to the formation of the 
jets simulated by Couch et al. (2009) in their vlml2 model. In Figure [T] we schematically 
summarize the proposed model. 

The duration of the jets launching phase from accretion disks around super-massive 
black holes in AGN is much longer than the dynamical time in the vicinity of the black hole. 




(1) 
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Fig. 1. — A schematic presentation of the proposed explosion model starting after bounce 
(not to scale). The four panels span a time of about one second, (a) As a result of the station- 
ary accretion shock instability (SASI; e.g. Blondin & Mezzacappa 2007; but see Nordhaus 
et al. 2010 for suppression of this instability in 3D simulations) or some other stochastic 
processes, segments of the post-shock accreted gas (inward to the stalled shock wave) pos- 
sess local angular momentum, (b) The accreted angular momentum changes stochastically 
in magnitude and direction. Namely, an intermittent accretion disk is formed or an accretion 
disk with a varying direction. The disk direction, and hence the direction of the jets, changes 
on a rate of ~ 10 rad s~^. (c) The jittering jets penetrate through the gas close to the 
center, but cannot penetrate beyond a typical distance of ~ 500 — 3000 km (eq. |9]). As 
the jet's axis changes its direction, no fresh jet's material is supplied to push the jet's head 
forward, as depicted by the short arrows detached from the center. The shocked jets form 
two hot bubbles. The bubbles accelerate the material around them in a wide angle (shown 
as short-wide arrows in the upper bubble only), forming slow massive wide (SMW) jets. 
This panel shows a larger volume than panel b. (d) The inflation of several pairs of bubbles 
eventually eject the entire region outside ~ 500 — 3000 km. This explodes the star. Unlike 
many other models, in the proposed model the stalled shock is not revived, but rather the 
gas within r < 500 — 2000 km, that amounts to ~ 0.1 — O.5M0, is accreted during the entire 
jet-launching phase. This panel shows a larger volume than panel c. 
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This is not the situation here. In our core collapse scenario the total accretion phase lasts 
for ~ 1 s, which is the typical time for the collapse of material from thousands of km. The 
intermittent accretion disk changes its axis direction on a time scale of ~ 0.1 s. The Keplerian 
period at the surface of the NS is ~ 0.003 s, such that the ratio of the typical time variation 
of the disk to the Keplerian period is only ~ 40. This implies that the accretion disk might 
not have time to completely relax and lose its internal energy and form a thin accretion disk. 
This non-relaxed behavior might further contribute to variations in the direction of the jets, 
and might increase the fraction of the accretion energy that is channeled to the outflow (jets) 
kinetic energy rather than to radiation. 

3. ANGULAR MOMENTUM CONSIDERATIONS 

In our model we assume that an intermittent accretion disc residing in the region 
20 km < r < 50 km is formed at the final stages of the collapse. There are two hmit- 
ing cases we now discuss. 

In the first limit, the total angular momentum of the accreted mass is very low, i.e., 
the average specific angular momentum is much lower than that of a Keplerian orbit about 
the radius of the NS. In that case, we assume, based on Ott et al. (2009), Blondin & 
Mezzacappa (2007), and Rantsiou et al. (2010) that the angular momentum is rapidly 
changing in magnitude and direction. As discussed above the value of the stochastically 
varying angular momentum amplitude can be large enough to form an accretion disk for a 
very short time of ~ 10 — 100 ms. The direction of the angular momentum stochastically 
changes as well. This stochastically variation of the angular momentum is not an independent 
'random walk' process. The reason is that the total angular momentum cannot be larger 
than that of the initial angular momentum of the core. We even expect the accreted angular 
momentum to be smaller than the initial angular momentum of the accreted mass. The 
reason is that the jets that are formed in each disk episode carry some angular momentum, 
and even a substantial fraction of the angular momentum in some models (e.g., Ferreira et 
al. 2006). The conclusion is that the NS will spin slower than what would be expected if 
the accreted mass would retain its angular momentum. We note that the small amount of 
accreted angular momentum in our model is along the results of Rantsiou et al. (2010) that 
found that the instability at shock does not lead to a rapidly spinning NS. We require the 
typical amplitude of the variations in the specific angular momentum to be large, but expect 
the total accreted angular momentum to be lower than the initial one. Our model, therefore, 
fits the finding that NS are born as slow rotators (e.g., Faucher-Giguere et al. 2006). 

The sum of the accreted angular momentum to almost zero, crudely requires the direc- 
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tion of the jets to rotate over 27r during the accretion process of ~ 1 s. This is compatible 
with, and one of the reasons for, taking ^ 10 rad s~^. 

Because of the stochastic accretion process the disk is not completely relaxed, and 
relative motion within the disk can increase the viscosity. For a disk viscosity parameter 
a = 1 the viscous timescale is (e.g., Dubus et al. 2001) 

. 5! . 10 f ^) " (^) " (^Y' (J^^) "'^ n,= (2) 

u V 0.2 y V 0.2 J VsOkmy \1AMqJ ^ ' 

where v is the viscosity of the disk, H is the thickness of the disk, Cs is the sound speed, 
is the Keplerian velocity, and Mns is the mass of the newly formed NS or black hole. As 
discussed in section |2] the intermittent disk is not expected to completely settled into a thin 
disk. This is the reason for taking H/Ra ~ 0.2 and Cg/v^ ~ 0.2 rather than a smaller value. 
We conclude that the viscosity is sufficient for accretion episodes that last ~ 10 — 100 ms. 

In the other limit the initial angular momentum of the accreted mass of the progenitor 
core is large, such that the disk maintains its general direction, and the direction of the 
jets is constant. This is probably the case in gamma ray bursts. If wide jets that are very 
efficient in expelling gas (Sternberg et al. 2007) are formed, such as in the magnetic model 
of Dessart et al. (2008), then small jittering will be sufficient to explode the star and leave a 
NS behind. Such a small jittering is expected from the shock instability. If narrow jets are 
formed, the small jittering will not be sufficient, unless more mass is accreted and launched 
into the jets. In such black hole will be formed. 



4. THE NON-PENETRATION CONDITION 

Let the fast jets from the inner disk zone have a mass outflow rate in both directions of 
M/, a velocity f/, and let the two jets cover a solid angle of 47r5 (on both sides of the disk 
together). The density of the outflow at radius r is 

_ 

From the numerical results for a 15Mq star at the beginning of the collapse as presented 
by Wilson et al. (1986) and Mikami et al. (2008), we find that the collapsing material density 
profile can be approximated by the power law 

Psir) = Por-^ = 1.3 X 10^° (^—L_y^' g cm'^, for 30 < r < 10^ km. (4) 
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The mass of the infaUing gas according to equation in the range 30 < r < 300 km, 
30 < r < 1000 km, and 30 < r < 2000 km, is 0.19Mq, O.35M0, and O.48M0, respectively. 
In our model, most of this mass is accreted to the newly formed NS star, bringing its mass 
from ~ I.IMq to a typical final mass of ~ 1.2 — I.GMq in most cases. This is the mass range 
where most NS studied by Kiziltan et al. (2010) reside. 

The head of each jet proceeds at a speed Vh given by the balance of the pressure exerted 
by the shocked jet's material with that of the shocked surrounding gas. Assuming a highly 
supersonic motion and neglecting the inflow velocity of the inflowing gas, this equality reads 
Psvl ~ Pf{vf — VhY, which can be solved for vj/vh 

^-i-f^V''=py^f^^V''. (5) 

\PfJ \ Mf j 

For typical parameters of our problem we find 

^ ^ 2 6 f ^ ( ^ ( ' \ + 1 (6) 

Vh ■ yO.O3M0 s-iy VO-01/ VlO^kms-V V 1000 km/ ' ^' 

where we took a typical mass of rjMg^cc — O.O3M0 to be blown into the two jets within a 
typical time of 1 second. We note that the accretion rate onto the newly formed neutron 
star can be quite high in the first ~ 0.3 s. In a ISM© it can be Mace — 0.3Mq s~^ after 0.2 s 
(J. Nordhaus, private communication 2011). The exponent —0.35 in the last term comes 
from 1 — 0.5a for a = 2.7. The time required for the jets to cross the surrounding gas at a 
distance r from the center is given by 

t ^ ^ ^ 04 ( ' ) ( ^ (^] ( ^ s (7) 

^ Vh ' V 1000 km/ yO.O2M0 s-M VO-01/ VlO^kms-i/ ' 

If the inflow velocity of the ambient gas is considered, the jets' head proceeds at a somewhat 
lower velocity than that given by equation (E]). On the other hand, the interaction of the 
jets with the ambient gas can further coUimate the jets, hence reducing somewhat the value 
of 6 and increasing the jets' head speed. 

As mentioned in section [H for the efficient operation of the penetrating jet feedback 
mechanism in expelling the ambient gas, the relative direction between the narrow jets and 
the ambient gas must be changing, such that the jets do not manage to drill a hole and 
escape with their energy. This implies that before a jet crosses a distance ~ r, its direction 
should be changed on an angle larger than its width. The time that a jet crosses its width 
Dj is Ts = Dj/{(j)r). The width of a jet at a distance r from its source is Dj = 2r sin a, where 
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OL is the half opening angle of the jet. For a narrow jet we have sin a ^ a ~ (25)^/^, and 



(/)r VO-01 / \ 10 s- 



-1 



.,= -i^ = 0.03(^) s. (8) 



From the demand for efficient energy deposition, jg < tp, we get the condition 

(\ -1.54 / . \ 0.77 

where we used assumption 4 of section [21 Mj = rjM^cci in equation ([9]). 

The main conclusion from this derivation is that the jets will stop rapid expansion (stop 
penetration) at a distance of ~ 10^ km (or ~ 500 — 3000). This is the region from where 
hot bubbles will be formed. They will accelerate the gas around them to form SMW jets, 
and explode the star. It is important to note that this occurs outside the stalled shock. In 
our model the stalled shock is not revived, but rather the material from its vicinity continued 
to be accreted and feed the accretion disk that launches the jets. We note that the radius of 
energy deposition does not depend on the half opening angle of the jets a, but does depend 
quite strongly on the rate at which the jets' axis direction is changing. 



5. THE ROLE OF NEUTRINO COOLING 

We compare the neutrino cooling time with the jets' expansion time. The dominant 
process is pair production e+ + e~ ^e,n,T + ^e,n,T- Based on the results of Schinder et al. 
(1987) and Itoh et al. (1989, 1996) we approximate the neutrino specific cooling rate by 

- 5 X 10^4 (^^^^^ ergcm-^s-^ (10) 

We check the cooling at the radius where the jets stop penetrating and start to infiate 
the bubbles. The jets stop penetrating at ~ 10^ km (eq. [9]). At this radius the jets' head 
proceeds at a speed of t;^ — 3 x 10^ km s~^ (eq. [6]). The typical time for the jets to cross the 
radius is ~ Vg/vh — 0.03 — 0.04 s. Each bubble is infiated by the shocked jets' material at 
a supersonic speed Ub > Cs — 10^ km s~^. We therefore take the bubbles infiation time to 
be tb ~ Tg/vh — 0.03 — 0.04 s. The volume of the bubbles infiated on both sides is 



Vb = 2x —(tbUbY ^ 2.3 X 10^'' —r^ r — — cm:'. (11] 

3 ^ ^ VIO^ km s-i/ VO.03 s ' ^ ' 
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The bubbles occupy ~ 0.1 of the volume when the jets reach ~ 10^. The total energy 
deposited into the bubbles is 



^50 ^ ?7 \ / Mace \ f Vf V ( h 



We emphasize that this is only the energy deposited in the first ~ 0.03 s. The jets are active 
for ~ 1 s. The energy in the bubbles is dominated by radiation. We can approximate the 
temperature by 

P \ 1/4 / + \ -1/2 



-1.5x10^" K, (13) 

^ aVj VO.03 s ,' ' ^ ' 



with the other parameters as used above. 



From equation ffTOj) we find that the total energy carried by neutrinos from the bubbles 
when the bubbles are inflated is 

/ 4. \ -0.5 

^--IO'MtTT^ tb>0.01. s. (14) 

\0.03 s / 

Bubbles that are inflated earlier lose more energy. For the above parameters, the ratio of 
energy carried by neutrinos to that deposited into the bubbles is 

r^^ ( tb ^ ^'^ 



Neutrino cooling is only important at very early times. The conclusion is that neutrino 
cooling is not a problem to the model. 



6. SUMMARY 

We examined a mechanism by which narrow jets launched by the newly formed neutron 
star (NS) or black hole (BH) at the center of a core collapse supernova (CCSN) lead to 
the supernova explosion of the star. Specifically, the mechanism converts fast narrow jets 
blown by the accreting newly formed NS, to slow massive wide (SMW) jets, similar to 
those simulated by, e.g.. Couch et al. (2009), as schematically shown in Fig. [H The wide 
outflow explosively expel the rest of the stellar mass. The basic assumption is that after 
a compact NS is formed, the final accretion stage results in the formation of an accretion 
disk that launches two fast narrow jets, similar in many respects to the way young stellar 
objects launch jets as they accrete mass during their final growth phase. The other basic 
assumptions are listed in section [2l The basic requirement for the jets to deposit their energy 
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into the infalling mass and turn the inflow to an outflow, is that the jets do not penetrate 
through the surrounding infalhng gas to too large distances. For CCSNe the jets should 
not penetrate beyond ~ 3000 km. If jets do penetrate, they carry most of their energy 
to distances of r > 10^ km, allowing too much mass to be accreted onto the NS, hence 
forming a BH. When the jets' axis does not change its direction and the surrounding gas 
has no transverse velocity, the jets drill their way out (Paper 1). For not penetrating the 
infalling gas the jets should encounter new material before they manage to drill their way 
out. Namely, the typical time for a jet's axis to cross the jet's width (eq. [H]), or for the 
infalling gas to cross the jet's width by transverse motion at radius (paper 1), should be 
shorter than the penetration time (eq. [7]) at radius r^: Tg < tp. In paper 1 the transverse 
motion of the inflowing gas was considered. Here we have studied the more realistic case 
where the jets' axis rapidly changes its direction: jittering jets. The condition < tp leads 
to equation 

This mechanism where the negative feedback is based on the jets' non-penetrating condi- 
tion is termed the penetrating jet feedback mechanism. It was suggested to operate in forming 
SMW jets (Soker 2008) and to regulate the growth of the super massive black hole during 
galaxy formation (Soker 2009; Soker & Meiron 2010). Namely, we suggest that NSs (or BHs) 
at the center of CCSNe shut off their own growth and expel the rest of the mass available 
for accretion by the same mechanism, via jets, that super massive black holes shut off their 
own growth, as well as that of their host bulge, in young galaxies. Another condition for the 
operation of the penetrating jet feedback mechanism is that the hot bubbles non-adiabatic 
cooling time is long. In the present case the non-adiabatic cooling is due to neutrinos. The 
neutrino cooling was studied in section |5] and was found to be very small. In most cases the 
jittering jets will lead to a spherical explosion. However, in some cases where the angular 
momentum of the accreted mass is relatively high, there will be a preferred direction for the 
varying jets' axis. In such cases the explosion will be bipolar. Such cases can explain the 
evidenced of asymmetry in some core collapse supernovae (e.g., Leonard et al. 2006; Wang 
et al. 2003). 

Our model has two signiflcant differences from most models for exploding CCSNe. (i) In 
our model neutrinos play no role in exploding the star. They can play a role in the accretion 
and launching of the jets, {ii) In our model the stalled shock at r ~ 300 km is not revived. 
To the contrary, it is required that material in the vicinity of the stalled shock is accreted 
and forms the intermittent accretion disk that launches the two jets. 

We note that, contrary to intuitive expectation, when the angular momentum of the 
collapsing core is very large the explosion is not efficient. This is because the jets' axis is 
constant, and the jets penetrate to large distances. Wide jets as in the simulations of 
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Dessart et al. (2008) can make the explosion more efficient. In most cases with large angular 
momentum we expect a large amount of mass to be accreted, and forms a BH. The two jets 
that are launched by the BH might form a Gamma Ray Burst (GRB). 

We thank Jason Nordhaus and the anonymous referee for helpful comments. This 
research was supported by the Asher Fund for Space Research at the Technion, and by the 
Israel Science Foundation. 
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